Gamma flashes from relativistic electron-positron plasma droplets 
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Ultra-intense lasers are expected to produce, in near future, relativistic electron-positron plasma 
droplets. Considering the local photon production rate in complete leading order in quantum elec- 
trodynamics (QED), we point out that these droplets are interesting sources of gamma ray flashes. 
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Fastly progressing high-intensity laser technology [l|, Q 
offers the perspective to create relativistic plasmas com- 
posed of electrons (e~), positrons (e"*") and photons (7) 
under laboratory conditions. Accordingly simulations 
have been performed, e.g., in Ql for special scenar- 
ios. Ideally would be the creation of an equilibrated 
e"'"e~7 plasma of sufficient size exceeding the mean free 
path of photons. Such a state of matter with properties 
reviewed in is governed by thermo-field theory ap- 
plied to QED and also in Q by kinetic theory. If such a 
plasma could be produced by table-top laser installations, 
it would represent the QED pendent to the QCD quark- 
gluon plasma currently investigated with large-scale ac- 
celerators 0] • Relativistic plasmas have some interesting 
specific features, such as collective plasmon or plasmino 
excitations 0] with accompanying van Hove singularities 
. Electron-positron plasmas are also interesting for as- 
trophysical scenario Temperatures of about T ~ 10 
MeV, in a laser-generated plasma, would open, further- 
more, channels for muon or pion production [ll|. While 
these channels are interesting for their own right (see 
[1^ for another avenue for laser driven particle produc- 
tion), a hot e"'"e~7 plasma droplet, exploding after the 
creation process by the enormous thermodynamic pres- 
sure, is a source of 7 radiation. Such short 7 fiashes may 
be of future use in ultra-fast spectroscopic investigations 
of various kind. 

In this note, we are going to present estimates of the 
photon spectrum emerging from a e^e~^ plasma droplet 
with temperatures in the 10 MeV range. In contrast 
to often employed particle-in-cell (PIC) simulations (see 
3 for a study of the start-up phase of a similarly hot 
plasma), we base our considerations on results from QED 
thermo-field theory. The expansion dynamics is treated 
in a schematic way, as our emphasis is on the photon 
emission characteristics. 

The calculation of the local spontaneous photon emis- 
sion rate from a QED (electron-positron) plasma has 
been outlined in [l3[ to complete leading order in elec- 



tromagnetic coupling a by including two-loop order. The 
result, at the heart of our note, may be presented as 
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where np (E-y) = {ex.p{E-y/T) + l)-i denotes the Fermi 
distribution, and the dynamically generated asymptotic 
mass squared of the electron is = 2tojj^. The thermal 
mass squared is given by m^j^ = with = Ana 

(The considerations apply for T > mg± which are dif- 
ferent from the non-relativistic case where the relevant 
scales are given by the masses of plasma particles and the 
temperature Q.) E^ = u ■ k is a short hand notation for 
the Lorentz scalar product of the medium's four-velocity 
u(t, x) and the photon's four-momentum k. This rate in- 
cludes 2 2 processes from one loop jl3l, viz., Compton 
scattering and pair annihilation which generate the lead- 
ing logarithmic contribution (first three terms in ([1])). In 
addition, the rate also includes the inelastic processes [HI 
from two-loop like bremsstrahlung (fourth term) and off- 
shell pair annihilation (fifth term) with the correct in- 
corporation [1^, [1^ of the Landau-Pomeranchuk-Migdal 
effect that limits the coherence length of the emitted ra- 
diation. A common feature of the latter two processes 
is that they have off-shell fermion next to the vertex 
where the photon is emitted, and the virtuality of the 
fermion becomes very small if the photon is emitted in 
forward direction. However, contrary to the one-loop di- 
agrams, the singularity is linear instead of logarithmic, 
and it brings a factor T^/m^j^ making it the same order 
in a as that of one-loop. The functions Cs are indepen- 
dent of a (because the relative importance of scatter- 
ing with a photon to an another charged particle in the 
plasma is essentially given by the ratio m^/m|), where 
niD is the Debye screening mass) but are nontrivial func- 
tions of y = Ej/T determined by a set of integral equa- 
tions |l3l [l6t. The numerical results for electron-positron 
plasma are described quite accurately by phcnomenolog- 
ical fits in the domain 0.1 < y < 30 which is sufficient for 
our consideration. The C22 is same in QED and QCD 
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because the interference effect cancels out and is given as 

C22(y) ^ OMly-^ - 0.3615 + LOle^i-^^^, (2) 

whereas the other two coefficients due to inelastic pro- 
cesses in QED are fitted together yielding 



Cb(y) + Ca(y) 



2aln(6+ 1/y) 

yC 
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where a = 0.374958, b = 0.431304, c = -0.05465, d = 
0.157472 and / = 1.73085. 

In Fig.[l] the scaled differential photon rate is displayed 
as a function of scaled photon energy. A comparison 
reveals that the low-energy photon rate is suppressed in 
the QED plasma relative to a QCD plasma (consisting of 
quarks, anti-quarks and gluons, all strongly interacting). 
The reason is that a photon has a non-negligible medium 
induced thermal mass, while in QCD its thermal mass is 
suppressed relative to the quarks by where g is 

the strong coupling in QCD. 
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FIG. 1: Differential photon emission rate (dashed) as well 
that of energy weighted (solid) for an equilibrated, symmetric 
e^e~7 plasma as a function of photon energy scaled with 
temperature. 



The damping of a photon is related to the absorption 
of it in the medium. The damping rate can be obtained 
from the photon production rate in ([T]) by the principle 
of detailed balance [13, 113, UM as 
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where the prcfactor, (27r)'^/4, is a matter of defini- 
tion [isj . The mean free path, A^, of a photon in a 
medium is inversely related to the photon damping rate 
as = ^1^1- In Fig- [2] the mean free paths of pho- 
tons arc displayed as a function of temperature of the 
electron-positron plasma for a wide range of photon en- 
ergies. This figure reveals very interesting features of the 
e"'"e~7 plasma droplet. Apparently at initial tempera- 
ture, it appears to be completely opaque as discussed in 
Refs. 0,|Tl|. However, it may not be so forever as we will 
see below. For a given photon energy the mean free 



E^=5 MeV 
E.^= 1 O MeV 
E,=25 MeV 
E,= 50 MeV 



T (MeV) 

FIG. 2: Mean free path of photons of energy E^, in an equi- 
librated, symmetric e^e~7 plasma as a function of tempera- 
ture. 



path increases with decreasing temperature implying a 
dilute system after some time. On the other hand, for a 
given T it increases with energy, i?^, making the system 
more transparent to the low energy photons depending 
upon the size of the plasma droplet. However, the real- 
istic scenario would be to consider the expansion of the 
droplet due to thermodynamic pressure. 

Wc consider the expansion with spherical geometry of 
radius R{t) = Rq + v{t)t, where Rq is the initial size of 
the droplet, v{t) is the surface velocity of the droplet 
assumed as v{t) = rjt, where 77 is a constant having 
dimension c^/ [Length]. We adopt a linear flow profile 
v{t,r) = v{t)r/R{t). The space-time integrated radia- 
tion spectrum follows from 

E^- — = 2 dt drr^ 
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where L(r , (j)) = [R^ (t) — sin^ (/))"'" — r cos (f) is the dis- 
tance [19[ the photon travels after it is created at a point 
r with an angle (f> relative to the radial direction. This 
expression takes into account that the radiation can only 
be emitted from such positions r and in such directions 
(j) where the distances to the surface R{t) is shorter than 
the mean free path A-y(<). The temperature at a time 
instant t is calculated by employing the energy conserva- 
tion relation: E{t) = Eq- £^™(t) - E^'^it), where Eq is 
the initial energy of the droplet obtained from the total 
initial energy density, e(To) = Utt^T^/GO ~ 3.8 x 10^^ 
Jm-3 at To = 10 MeV. We note that Eq (oc T(fi?i] ) de- 
pends crucially on the initial droplet size. Eq ~ 13 kJ for 
a droplet of radius Rq ~ 2 ??m with the same reasoning 
as in Ref. [ll[ . (It is an interesting and challenging topic 
how the laser pulse energy can be deposited in a smallest 
possible spatial and temporal volume. The scenario in 
Ref. d, [U envisages, e.g., an initial stage of compressed 
foil with (1 - 2) nm thickness and Tq ~ 10 MeV). Now, 
at any instant of time t, one can obtain the energy of the 
emitted photon as E^'^{t) = J dE^ (E^dN/dE^) from 
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([5]); and the kinetic energy of the droplet as E^"^{t) = 
4tt J dr e{T{t))jv{t,r), with 7 = l/{l-v^{t,r)/c^Y^^. 
The evolution is truncated at time tf at which T cools 
down to 1 MeV. For lower temperatures, the electron 
rest mass becomes a relevant scale which modifies the 
employed rate. As we are interested in the hard part of 
the emission spectrum, > 1 MeV, the late evolution 
with T < 1 MeV is not essential for our purposes. 




10 15 20 25 

t (n m/c) 



FIG. 3: Temperature as a function of time for rj = 0.03 c^/nm 
(solid) and 0.01 c^/nm (dashed). 



In Fig. [21 the evolution of T{t) is displayed for two 
surface velocities. The expansion lasts up to ~ 
(20 — 30) nm/c. A realistic three-dimensional expansion 
dynamics can be obtained from relativistic hydrodynam- 
ics from given initial conditions, in particular phase space 
distribution. This deserves separate dedicated investiga- 
tions which we postpone. A relevant scale is the velocity 
of sound (l/v^ in an ideal relativistic plasma) at which 
a rarefaction wave travels from the surface towards the 
center is compatible with our consideration. We note 
that the evolution depends strongly on the system size 
but weakly on the surface expansion velocity. 
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FIG. 4: Mean free path of a photon with E-y — 25 MeV as a 
function of time for 77 = 0.03 c^/nm (solid) and 0.01 c^/nm 
(dashed) . 



In Fig. m we display the evolution of A.^ (i) with same 
initial conditions as above. The mean free path of the 
photon for the expanding droplet depends on the velocity 
profile. This in turn determines the region from which the 



photons may escape. The system becomes more dilute 
and transparent in later time. Inspection of Fig. [5] also 
unravels that harder photons have a larger optical depth, 
as mentioned above. 
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FIG. 5: Evolution of £(t) (sohd), _B°'"(t) (dashed) and £'''"(*) 
(dash-dotted) for r] — 0.03 c^/nm and 0.01 c^/nm. A larger 
value of tf corresponds to a smaller value of 77. 



In Fig. O we exhibit the evolution of various energies 
with same initial conditions as above. Initially, E^^'^ in- 
creases with the increase of time but then decreases as 
the size of the system becomes larger and the local en- 
ergy density drops. The integrated energy of emitted 
photons increases with time indicating that the most of 
the system's energy is attributed to photon production. 
Thus, the total energy remaining in the droplet decreases 
with time. The surprising fact is that cooling by photon 
emission is very efficient, i.e. the initial energy of the 
plasma droplet is converted essentially into hard photon 
radiation. This limits the time span the droplet spends 
in a hot stage. In addition, e"*" and e~ leakage reduces 
the life time of the hot era further. The droplet, after 
equilibration, appears a gamma flash. Of course, there 
are relaxation processes, like heat conductivity etc. The 
present estimate assumes that these are short compared 
with the expansion time scale. 

In Fig. [SI we display the energy weighted photon spec- 
tra. As seen there is an enormous amount of photon pro- 
duction from the electron-positron plasma droplet. Note 
that the time integrated photon spectrum scales roughly 
with the initial volume or energy Eq, respectively. For 
instance, a one-order change in initial radius would cor- 
respond to a three-orders of magnitude change in photon 
production. 

We now briefly discuss an asymmetric plasma with 
non-zero net electric charge. The chemical potential 
fi steers such an asymmetry. The ratio of electron to 
positron densities is approximately oc exp{2^/r}, while 
the net density behaves like oc ^T^ in leading order. 
That means, a noticeable difference in and e~ den- 
sities translates into fi > T. fj. enters the photon rate in 
(dD through = e'^T^il + fi'^/{TT^T^))/4:, which takes 
into account the number of charge particles available for 
various processes. Also various Cs depend on fi through 
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FIG. 6: Energy weighted photon spectrum as a function 
of photon energy for surface expansion velocities with i] — 
0.03 c'^/nm (sohd) and 0.01 c^/nm (dashed). Radiation from 
the ignition stage is not included. 



the ratio ml^/ml - (1 + ny{TT^T^))/(l + i^i^ / {n^T'^)) 
which accounts for, as discussed earlier, the relative im- 
portance of scattering with a photon to a charge in the 
plasma. Therefore, it could only be important at fi > T. 

The present considerations are armed with a state-of- 
the art photon cmissivity. The dynamics, in contrast, is 
rather schematic: constant temperature over the fireball, 
linear velocity profile, prescribed surface velocity and ve- 
locity profile etc. Envisaged calculations within relativis- 
tic hydrodynamics including gradients and details of ra- 
diation transport will make the scheme more realistic. 



Clearly, chemical undersaturation [3|, |20[ (which would 
correspond to a decrease of the initial energy density for 
the same initial temperature) parameterized by fugaci- 
ties need to be included along the lines of [2l[ to arrive 
at a detailed picture. The effect of impurities, like ions or 
the original material where the plasma is formed from, 
should be considered too, as other expansion patterns, 
e.g, for the geometry envisaged in Rcf. 3]. 

In summary we point out that a hot e~^e~j plasma 
droplet created by ultra-intense and short laser pulses 
is a technologically interesting source of 7 flashes with 
continuous spectrum (and not only a flash of 511 keV 
radiation as argued in Ref. 3) ranging up to energies 
being a multitude of the initial temperature. To have an 
estimate of the efficiency for transforming the incoming 
eV laser photons into hard gamma radiation, we mention 
that the integrated power in the 7 flash is comparable to 
the total energy of the initial plasma droplet. Hereby 
we considered only photons with > 1 MeV. In the 
considered scenario, the other energy, after disintegra- 
tion, resides in collective (radial) flow of electrons and 
positrons and soft photons. 
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